1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
| library(MASS)
library(class)
x0<- MASS::mvrnorm(n=10, mu = c(1,0), Sigma = diag(2))
x1<- MASS::mvrnorm(n=10, mu = c(0,1), Sigma = diag(2))
make.data<-function(x0,x1,n,m,Sigma1 = diag(2)/5,Sigma2 = diag(2)/5){
r<-length(x0[,1])
x<-matrix(0,n+m,2)
for(i in 1:n){
s<-sample(1:r,1)
x[i,]=MASS::mvrnorm(n=1, mu = x0[s,], Sigma = Sigma1)
}
r<-length(x1[,1])
for(i in 1:m){
s<-sample(1:r,1)
x[i+n,]=MASS::mvrnorm(n=1, mu = x1[s,], Sigma = Sigma2)
}
y<-matrix(c(rep(0,n),rep(1,m)),n+m,1)
data<-data.frame(y,x)
return(data)
}
train.data<-make.data(x0,x1,100,100,diag(2)/5,diag(2)/5)
|
